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Abstract 

Particle production processes in the expanding universe are described 
within a simple kinetic model. The equilibrium conditions for a Maxwell- 
Boltzmann gas with variable particle number are investigated. We find 
that radiation and nonrelativistic matter may be in equilibrium at the 
same temperature provided the matter particles are created at a rate that 
is half the expansion rate. Using the fact that the creation of particles is 
dynamically equivalent to a nonvanishing bulk pressure we calculate the 
backreaction of this process on the cosmological dynamics. It turns out 
that the 'adiabatic' creation of massive particles with an equilibrium distri- 
bution for the latter necessarily implies power-law inflation. Exponential 
inflation in this context is shown to become inconsistent with the second 
law of thermodynamics after a time interval of the order of the Hubble 
time. 

PACS numbers: 98.80.Hw, 05.20.Dd, 04.40.Nr, 05.70.Ln 



1 Introduction 

Particle production processes in the early universe are supposed to have in- 
fluenced the cosmological history considerably. These processes are quantum 
in nature [|l|. Their backreaction on the cosmological dynamics, however, are 
frequently studied phenomenologically. As was observed by Zel'dovich and 
Murphy Q and lateron confirmed by Hu [Q] a nonvanishing particle production 
rate is equivalent to a bulk viscous pressure in fluid cosmology. This equiv- 
alence has been used recently in a series of papers [^, 0, ||8|, |^, [10|, 



pd[| , [^] , , , discussing different aspects of this effective viscous pressure 
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approach. Taking into account that the fluid dynamics of a gas may be de- 
rived from Boltzmann's equation for the one-particle distribution function, the 
problem has been addressed whether the effective viscous pressure approach is 
compatible with the kinetic theory of a gas with varying particle number |15|. A 
modified Boltzmann equation was proposed in which an additional source term 
describes the change of the one-particle distribution function due to particle 
number nonconserving processes, supposedly of quantum origin. For a sim- 
ple rate approximation of this source term and under the assumption that the 
characteristic time scale for interactions between the particles is much smaller 
than the time scale on which the number of particles changes, it turned out 
that the effective viscous pressure approach is compatible with kinetic theory 
in homogeneous spacetimes but not in inhomogeneous ones [15|. 

A modification of Boltzmann's equation implies modified conditions for col- 
lisional equilibrium. In the standard relativistic kinetic theory, i.e., without the 
mentioned additional source term, it is a well established result |17| that a 
simple gas can only be in 'global' equilibrium in spacetimes that admit a time- 
like Killing- vector, i.e., in stationary spacetimes. As an immediate consequence, 
it is impossible to characterize such a gas by a 'global' equilibrium distribution 
function in an expanding universe. The only exception is a gas of massless 
(m = 0) particles (radiation), where the corresponding equilibrium condition is 
less restrictive and requires only the existence of a conformal Killing-vector. 

As was shown in [15|, the conditions for the one-particle distribution func- 
tion to have a structure characteristic for collisional equilibrium are weaker in 
the case of a modified Boltzmann equation, i.e., a Boltzmann equation with an 
additional source term, than in the standard case. More specifically, for the 
mentioned effective rate model of the source term the one-particle distribution 
function takes its collisional equilibrium structure if the corresponding space- 
time admits a conformal Killing-vector for massive particles as well. A gas of 
particles with m > in the expanding universe may also be characterized by a 
distribution function with equilibrium structure, provided the number of par- 
ticles is not fixed but changes according to a specific rate, depending on the 
detailed equations of state. 

The present paper provides a detailed discussion of the conditions for col- 
lisional equilibrium of a simple gas with nonconserved particle number. Our 
main objective is to consider the possible implications of these solutions of the 
modified Boltzmann equation for the cosmological dynamics. 

Of course, the applicability of a kinetic approach to the early stages of 
the cosmological evolution is not straightforward. A gas, however, is the only 
system for which the correspondence between microscopic variables, governed 
by a distribution function, and phenomenological fluid quantities is sufficiently 
well understood. All the considerations of this paper refer to a model universe 
for which the kinetic approach is assumed to be valid. We hope that this 
idealized model nevertheless shares some basic features with our real universe. 

Under these prepositions we shall establish apparently surprizing links be- 
tween the equilibrium properties of a relativistic gas with nonconserved particle 
number and the spacetime geometry of the universe. The most striking feature 
is the circumstance that the collisional equilibrium of a classical Boltzmann gas 
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with varying particle number in some cases requires a power-law inflationary 
universe. 

The conditions for collisional equilibrium do not fully fix the particle pro- 
duction rate. The remaining freedom may be used to impose additional physical 
requirements. We shall discuss here two different cases and explore the impli- 
cations of each of them for the cosmological evolution. 

The first case will be that of 'adiabatic' particle production, characterized by 
a constant entropy per particle. The particles areassumed to be created with 
a fixed entropy. This additional assumption entirely determines the creation 
rate that is neccessary for given equations of state to keep the particles at 
equilibrium. This rate is the higher the more massive the particles are. For 
nonrelativistic particles it is just half the expansion rate. 

Now, a given creation rate is equivalent to a given effective viscous pressure. 
With the help of the latter it is possible to calculate the backreaction of the 
production process on the cosmological dynamics. For nonrelativistic particles 
this backreaction leads to power-law inflation with a behaviour of the scale 
factor R like R oc t^/^ 

in a homogeneous, isotropic and spatially flat universe. 
Despite of the equilibrium structure of the distribution function the generation 
of particles is accompanied by a nonvanishing entropy production. For massive 
particles the comoving entropy grows as t'^. 

The second case is connected with exponential inflation. There have been 
suggestions in the literature that a sufficiently high production rate of particles 
or strings might give rise to a de Sitter phase (|jl^, @, [^, [10|, [12|). We 
shall consider here the question whether this suggestion is compatible with the 
kinetic theory of a gas with variable particle number. To this end we use the 
above mentioned freedom in the production rate to impose the condition of a 
constant Hubble rate and clarify the consequences of this assumption. We show 
that exponential inflation becomes inconsistent with a non-negative entropy 
production, as required by the second law of thermodynamics, after a time of 
the order of the Hubble time. 

The paper is organized as follows. In section 2 we recall the previously de- 
veloped kinetic theory for particle production and discuss the general conditions 
for collisional equilibrium. Section 3 is devoted to 'adiabatic' particle produc- 
tion and its relation to power-law inflation. In section 4 we check an inflationary 
scenario where a de Sitter phase is generated by 'nonadiabatic' production of 
gas particles with a Maxwell-Boltzmann equilibrium distribution function and 
show that it violates the second law of thermodynamics after about one Hubble 
time. The final section 5 summarizes the main results of the paper. 
Units have been chosen so that c = ks = h = 1. 



2 Kinetic theory for particle production 
2.1 General kinetic theory 

According to the discussion in ||l^ we assume that a change in the number of 
particles of a relativistic gas should manifest itself in a source term H on the 
level of kinetic theory. 
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The corresponding one-particle distribution function f = f {x,p) of a rela- 
tivistic gas with varying particle number is supposed to obey the equation 



L [f] = p'fn -niP'p'^ = C[f]+H {x,p) , (1) 

where f {x,p) p^nf^dT,dP is the number of particles whose world lines inter- 
sect the hypersurface element n^dTi around x, having 4-momenta in the range 
{p,p + dp); i, k, I ... = 0, 1, 2, 3. 

dP = A{p)5 {p^Pi + m?) dP4 is the volume element on the mass shell p'^pi = —m? 
in the momentum space. A{p) = 2, if is future directed and A[p) = other- 
wise; dPi = —gdp^dp^dp^dp'^ . 

C[f] is the Boltzmann collision term. Its specific structure discussed e.g. by 
Ehlers will not be relevant for our considerations. Following Israel and 
Stewart we shall only require that (i) C is a local function of the distri- 
bution function, i.e., independent of derivatives of /, (ii) C is consistent with 
conservation of 4-momentum and number of particles, and (iii) C yields a non- 
negative expression for the entropy production and does not vanish unless / 
has the form of a local equilibrium distribution (see (|^). 

The term H{x,p) on the r.h.s. of (|l]) takes into account the fact that the 
number of particles whose world lines intersect a given hypersurface element 
within a certain range of momenta may additionally change due to creation or 
decay processes, supposedly of quantum origin. On the level of classical kinetic 
theory we shall regard this term as a given input quantity. Later we shall give 
an example for the possible functional structure of H{x,p). 

By the splitting of the r.h.s of eq.(|l|) into C and H{x,p) we have separated 
the collisional from the creation (decay) events. In this setting collisions are 
not accompanied by creation or annihilation processes. In other words, once 
created, the interactions between the particles are both energy-momentum and 
number preserving. For a vanishing H{x,p) eq.(||) reduces to the familiar Boltz- 
mann equation (see, e.g., |16], [18|, [l^, pO|). 

The particle number flow 4-vector A^* and the energy momentum tensor T'^^ 
are defined in a standard way (see, e.g., [^]) as 

N' = I dPp'fix,p) , = I dPp'p^f{x,p) . (2) 

The integrals in (Q) and throughout the paper are integrals over the entire mass 
shell p^pi = —m?. The entropy flow vector 5" is given by |17| 



5'^ = - |p"[/ln/-/]dP, (3) 

where we have restricted ourselves to the case of classical Maxwell-Boltzmann 
particles. 

Using the general relationship [18| 

y p'^i . . . .p''"//dP = y P"' • • -P"" L [/] dP (4) 
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and eq.(|T|) we find 

N;^a = J{C[f] + H)dP, T-;i = j p'^{C[f]+H)<iP , (5) 

and 

S% = - j\nf{C[f]+H)dP . (6) 

In coUisional equilibrium, which we shah assume from now on, In/ in (^) is 
a Unear combination of the colhsion invariants 1 and p". The corresponding 
equihbrium distribution function becomes (see, e.g., |]l6|) 



/°(x,p) = exp[a + M , (7) 

where a = a (x) and fia is a timehke vector that depends on x only. 
Inserting the equilibrium function into eq.(0) one gets 

P>a + /?(a;fe)Py]/° = ^(x,p) . (8) 

It is well known ([^], fl^) that for H(x,p) = this equation, which char- 
acterizes the 'global equilibrium', admits solutions only for very special cases 
in which a = const and is a timelike Killing-vector. It will be the main 
objective of this paper to clarify the consequences of ^ for a specific choice of 
the source term H{x,p) to be discussed below. 

With (0), the balances @ reduce to 



J. y 



J HdP , Tfk = J p^HdP . (9) 



In collisional equilibrium there is entropy production only due to the source 
term H. From eq.(^) we obtain 

S% = - J Hix,p)lnfdP , (10) 

implying S% = -aN% - paT%. 

With / replaced by f in (|) and (|), N", T"^ and may be split with 
respect to the unique 4-velocity u"" according to 

iV^ = nu" , T"^ = pu'^v!' + ph""^ , S"" = nsu" , (11) 

where /i'^^ is the spatial projection tensor h""^ = g""^ + u°'u^, n is the particle 
number density, p is the energy density, p is the equilibrium pressure and s is 
the entropy per particle. The exact integral expressions for n, p, p and s are 
given by the formulae (177) - (180) in [16|. 
Using (pT|) and defining 



T = ^jH{x,p)dP, (12) 

the first eq.(|^ becomes 



h + On = nT . (13) 
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It is obvious that T is the particle production rate. Similarly, with the decom- 
position ([TT| ) and the abbreviation 

e = - J p''H{x,p)dP , (14) 

the energy balance, following from the second equation (^, may be written as 

p+Qip + p)-Uat'' = , (15) 

where = is the fluid expansion. 
Introducing the quantity 

vr^-- 1-, (16) 



it is possible to rewrite the energy balance (15) as 



p + G(/) + p + 7r) = 0, (17) 

i.e., vr enters the energy balance in the same way as a bulk viscous pressure 
does. A rewriting like this corresponds to the introduction of an effectively 
conserved energy-momentum tensor 

= pu'u^ ^{p^T^)W^ , (18) 

instead of the nonconserved quantity T^^ of (0). This interpretation of the 



source term ( |14|) as an effective bulk pressure was shown to be consistent for 



homogeneos spacetimes ]15| . 

2.2 Effective rate approximation 

The quantity H (x,p) is an input quantity on the level of classical kinetic the- 
ory. It supposedly represents the net effect of certain quantum processes with 
variable particle numbers (see, e.g., [|l|) at the interface to the classical (non- 
quantum) level of description. Lacking a better understanding of these processes 



it was assumed in [15| that the influence of these processes on the distribution 
function f {x,p) may be approximately described by a linear coupling to the 
latter: 

H{x,p) = (:{x,p)f{x,p) . (19) 
Moreover, C was supposed to depend on the momenta only linearly: 

C = -^ + u{x). (20) 

C, or equivalently and r, characterize the rate of change of the distribution 
function due to the underlying processes with variable particle numbers. The 
restriction to a linear dependence of C on the momentum is equivalent to the 
requirement that these processes couple to the particle number flow vector and 
to the energy momentum tensor in the balances for these quantities only, but 



6 



not to higher moments of the distribution function. This 'effective rate approxi- 
mation' is modeled after the relaxation time approximations for the Boltzmann 
collision term (p^, |23], 123]). While the physical situations in both approaches 



are very different, their common feature is the simplified description of nonequi- 
librium phenomena by a linear equation for the distribution function in terms of 
some effective functions of space and time that characterize the relevant scales 
of the process under consideration. In the relaxation time approximation this 
process is determined by the rate at which the system relaxes to an equilibrium 
state. In the present case the corresponding quantity is the rate by which the 
number of particles changes. 

With @ and dH) for H{x,p) the source terms nV and in (|l3|) and (||) 
are given by 

7/ 71 

nT = -—N'' + u{x)M{x) = - + uM, (21) 
r r 

where M is the zeroth moment of the distribution function, M = J dPf {x,p), 
and 

- Uaf" = + „ (x) n.iV" = -P-un, (22) 

r [X) T 

respectively. The expression (|2^) for ensures that the source terms (pT|) 
and (^) depend on M, iV* and T*^ only, but not on higher moments of the 
distribution function. 

Using the specific structure (|l^) and ( pO|) of the source term H{x,p) in (^) 
the condition (|8|) becomes 

P"«,a + /3(a;fe)Py = f + ^' (23) 

where E = —Uap"'- Decomposing p"" according to p"" = Eu"" + Ae", where e" 
is a unit spatial vector, i.e., e"ea = 1, u°'ea = 0, the mass shell condition 
p"'Va = —Tn?' is equivalent to \^ = E"^ — m? . The conditions on a turn out to 
be 

d = - , hla^a = 0, (24) 
r 

while (5a obeys the equation for a conformal Killing vector 

I3{a;b) = ^(x)5afe , {3 i^a;b)U'' v!' = V , (25) 

implying ^ = —vjrn}. Different from the case H{x,p) = where a has to be 
constant in space and time, a is only spatially constant in the present case but 
changes along the fluid flow lines. From (^5[) it follows that z/ = for m = 0. 
The condition u = Q for m > 0, however, is equivalent to [3{a;b) = 0) the 
corresponding spacetime is stationary. 

Using 13a = Ua/T in yields f^, [|l|] 



T 1 

^ = — (26) 
T 3 ^ ' 

for the temperature behaviour in an expanding universe. It is remarkable that 
this relation which is a consequence of the conformal Killing-vector property 
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for /?" holds both for m = with u = and for m > in the case z/ 7^ 0. The 
validity of the same equilibrium relation (p6| ) both for radiation with H{x,p) = 
and for massive particles with H{x,p) 7^ has the interesting implication that 
radiation and matter in the expanding Universe may be in equiUbrium at the 
same temperature, provided the particle number of the matter component is 
allowed to change. It is well known, that for conserved particle numbers, i.e., 
H{x,p) = u = = 0, an equilibrium between both components is impossible 

As a consequence of the fact that Ua/T is a conformal, timelike Killing- 
vector, the 4-acceleration u"" may be expressed in terms of the spatial temper- 
ature gradient according to 

Ua = -(lnT),fe/it (27) 

Within the linear theory of irreversible processes this coincides with the condi- 
tion for the heat flow to vanish (see, e.g. [^6|). For a comoving observer in the 
rotation free case (p7|) reduces to[]T5| 



(T^/^),. = , (1^ = 1,2,3) . (28) 

The latter formula which states that the quantity T^J—qqq has to be spatially 
constant if particles with m > are produced, replaces Tolman's relation 
{T^f--gm),n = 0, according to which T^J—goQ is constant both in space and 
time in case Ua/T is a Killing-vector. 



The conformal Killing-equation (25) provided us with the temperature law 



(Pq). On the other hand, the behaviour of the temperature is determined ther- 
modynamically. The fluid equations of state may generally be written as 

p = p{n,T) , p = p{n,T) . (29) 

Differentiating the latter relation and using the balances ([T3|) and ( |T7|) one finds 

'dp/dT p + P 



T 



dp/dT vr 
dp/dT ^ Tdp/dT 



+ r 



dp/dT Tdp/dT 



(30) 



where dp/dT = {dp/dT)^ and dp/dT = (dp/dT)^. Relation (||) was first 
derived in 

In the specific case of a classical gas the equations of state are p = nT with 

i2r 



m 



ATTrri^T 

and p = ne{T) with 

Ki (^) , 



^ = ,^ x3 ( 7F ) exp [a] , (31) 



The quantities are modified Bessel- functions of the second kind Em- 
ploying the well-known differential and recurrence relations for the latter, one 
gets 

dp _ dp/dT 2 h fhy 
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where z = m/T and 



h = e + 



P 



m- 



(34) 



n K2{z) ' 
is the enthalpy per particle. 

By virtue of (16), (^) and (^), the temperature law ( |30|) may be expressed 
in terms of u and r. Because of dp/dT = p/T the r-terms cancel and we obtain 



T 
T 



vn 



dp/dT Tdp/dT 
Combining the latter expression with (p6|) yields 



n n 



T- 



1 



Idp 
3 dp 



1 



pM_ 

n n 



e 



(35) 



(36) 



The i^-part of the source terms (^) , (^Tj) and ( p2[ ) is fixed by the thermodynamic 
functions of the gas. The zeroth moment M of the classical distribution function 
is given by 

f . 



M 



exp [a\ 



(37) 



There exists an alternative formula for the function v that may be obtained 
from the Gibbs-Duhem equation 

dp = nsdT + ndjj, . (38) 

With [||] 

p + p n/i 

ns = — — — 

T T 

one finds 



p={p + p)^ + nT{^ 



(39) 
(40) 



Using here p = nT yields 



Applying (|13D with (|2l|) for ii/n, eq.(|2 
a = {ji/T) , wearrive at 



n 
n 



Pi 
p T 



(41) 



for T/T and the first of eqs.(24) for 



-d 

M 



3p 



3 



\-z^ 

K2 



e 



(42) 



The equivalence between the expressions ( pq ) and (^2|) becomes manifest if one 
uses the recurrence relations for the functions Kn (see, e.g., H^) to replace, 
e.g., Ki {z) in M in (|6|) by K2 and K^^. With v given either by (|3|) or (||), 
the conformal factor ^ in ( p5| ) is fixed as well: 



J_K2 

3m Ki 



K3 
''K2 



G 



(43) 



i.e., it is completely determined by the thermodynamic functions of the gas. 
In the limit 2: 3> 1 we have v ~ —mzQ/'i and a corresponding expression for 
the conformal factor. While v is uniquely determined by the condition (p^), 
the function r is still arbitrary. We are free to impose a further condition to fix 
r. 
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3 Adiabatic particle production and power-law in- 
flation 



3.1 The creation rate 



In H] as well as in earlier papers (§, |], @, Jill, 11) ^he particle 

production was assumed to be 'adiabatic'. This implies that all particles are 
amenable to a perfect fluid description immediately after their creation. The 
particles are produced with a fixed entropy, i.e., the entropy per particle s does 
not change. The corresponding condition s = determines r. 
Prom the Gibbs equation 



Tds = d- +pd- 
n n 

together with (^) and (15) we find 

nTs = Uat" - {p + p)T . 
According to (^5|) the condition s = is generally equivalent to 

Uaf" = {p + p)T , 



(44) 



(45) 



(46) 



relating the source term in the energy balance to that in the particle number 
balance. Using ( pl[ ) and (22) with ( |36[ ) in ( ^6| ) provides us with 



ldp\ i_ p±pM 



3 dp 



pM 
n n 



e 



Because of ( ^T] ) the particle production rate T turns out to be 



3 dp. 



e 



in this case, equivalent to 



TT 



3 dp 



{p + p) , 



(47) 



(48) 



(49) 



where we have used (16). In fact, expression (^) for the adiabatic particle pro- 
duction rate may be obtained from (|2^), (^), (4^) and (16) without explicitly 
knowing v and r. A separate knowledge of i' and r, however, will be needed in 
section 4 below. 

Inserting into (E^) the expression (|33[) for dp/ dp we have 



1 + 5- 



T 



e , 



(50) 



or an equivalent expression for vr. In the limiting case m = (radiation) with 
h = 4T, the relation ( ^0|) yields F = in agreement with the result of ||l5|] that 
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the adiabatic production of massless particles is forbidden if P'^ is a conformal 
Killing-vector. In the opposite limiting case z ^ 1 the enthalpy per particle is 



h 



Z+- + 
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-1 



yielding 



-e 



(51) 



(52) 



i.e., the production rate for massive particles in equilibrium is half the expansion 
rate. Equivalently, since the Hubble factor of the Robertson- Walker metric is 
H = 0/3, the characteristic time-scale for the production of very massive 
particles, neccessary for an equilibrium distribution of the latter, is 2/3 of the 
Hubble time in the case s = 0. 



3.2 Backreaction on the cosmological dynamics 

The main advantage of a fluid approach to particle production is the possibility 
to calculate the backreaction of this process on the cosmological dynamics. For 
homogeneous spacetimes this problem is equivalent to studying the dynamics of 



a bulk viscous fluid universe |15]. While generally the determination of the bulk 
pressure is an involved problem on its own we are in a much better situation 
in the present case. The effective bulk pressure here is completely fixed by the 
conformal Killing- vector conditions (p5[), the adiabaticity condition s = and 
the equations of state. Using (^) in the energy balance (p!7D, the latter takes 
the form 

p + pdp 

with dp/ dp given by (p3|). While for m = 0, equivalent to p = p/3, with 
TT = we recover the familiar behaviour for radiation, the opposite limiting 
case, z 3> 1, yields 

P = -\hp, (54) 

or, with H = R/R, where R is the scale factor of the Robertson- Walker metric, 

p oc R-^/^ . (55) 

The energy density decreases less than without particle production {p oc i?""^). 
A parallel statement holds for n. From (|l^) with ( ^2|) we find n oc R^^^^ instead 
of n oc R-^ for F = 0. 

We recall that the temperature behaviour is given by (p6|), irrespective of 
the equations of state. 

Restricting ourselves to a homogeneous, isotropic and spatially flat universe 
with 

3-^ = Kp , (56) 
where k is Einstein's gravitational constant, and 

H-2p' ^^^^ 
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where the latter follows by virtue of 3R/R = — {k/2) {p + 3p + Sir), we find 
with the help of (|5^ ) that the scale factor R behaves like 

R oc t^/^ , (58) 

instead of the familiar R oc t^/^ for p oc R~^, i.e., for F = 0. 

If massive particles are dynamically dominating, the adiabatic production 
rate for these particles that is neccessary to keep them governed by an equi- 
librium distribution function backreacts on the dynamics of the Universe in a 
way that implies power-law inflation. In other words, massive particles with 
nonconserved particle number and fixed entropy per particle are allowed to be 
in collisional equilibrium only in a power-law inflationary universe. 

The backreaction is largest for z ^ 1 and it vanishes for m = 0. In the 
latter case the familiar R oc t^/^ behaviour for radiation is recovered. There 
exists an intermediate equation of state yielding R (xt with R = 0. The obvious 
condition for i? = is 

p + 3p + 37r = . (59) 
Using the expression (^) for tt, this condition reduces to 

1^ = ^. (60) 

op p + p 

Inserting here (|33t), we obtain the following cubic equation for h/T: 

^)'-5(^)'-(.^-3)A_2 = 0. (61) 

Together with (|3j) the latter equation determines the criticalvalue Zcr that 
characterizes the case R = 0. The numerical result is Zcr ~ 9.55. We have 
i? > for z > Zcr while z < z^r corresponds to i? < 0. 

In terms of a '7'-law, i.e., p = (7 — l)/9, the corresponding critical 7 value 

7cr is 

- - ~ 1.09 . (62) 



■ K2{zcr) 

There is power-law inflation for 1 < 7 < 7cr. The case 7 = 7cr corresponds to 
R = Q. 



3.3 Entropy production 

It is a specific feature of our approach that a gas, kept in equilibrium in the 
expanding universe by a nonvanishing particle creation rate exhibits a nonvan- 
ishing entropy production density 

S% = nVs , (63) 

with r from ([48|). While there is no entropy production due to conventional 
dissipative processes, it is the particle production rate itself that is connected 
with an increase in the entropy. 
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The comoving entropy S is defined by S = nsR^. This quantity may vary 
either by a change in the entropy per particle s, or by a varying number = 
nR^ of particles. The first possibility was excluded here by the requirement 
(^).The second one, an increase in the number of m > particles, however, 
was necessary for the particles to be governed by an equilibrium distribution 
function. Consequently, this kind of equilibrium even implies an increase in the 
comoving entropy. 

While the latter statement sounds unfamiliar and does never hold indeed 
for systems with conserved particle numbers, it is a natural outcome for a 
nonvanishing creation rate T. 

A further unfamiliarity, related to the previous one, is that one has nonzero 
entropy production in the corresponding power-law inflationary phase. In the 
standard inflationary scenarios the phases of accelerated expansion are not ac- 
companied by an increase in E and all the entropy is produced during a subse- 
quent reheating period. In our approach one finds from (^), ( 
the limiting case of very massive particles z ^ 1, that S ~ t^, i.e., the comoving 
entropy grows quadratically with the cosmic time. 



52) and (||)in 



4 Exponential inflation 

The functions ly and r characterizing the particle productionprocess were con- 
sidered to be input quantities on the level of classical kinetic theory. Only the 
specific expression (^6|) for however, guarantees the fulfilment of the equilib- 
rium conditions - ([25|). Once the remaining freedom in r is fixed by some 
physical requirement like that of the 'adiabaticity' of the creation process in the 
preceeding section, the dynamics of our model universe is uniquely determined. 
While the 'adiabaticity' condition appears 'natural' it is not the only possible 
choice. In the present section we use the freedom in r to impose the condition of 
a constant Hubble rate, i.e., the condition for exponential inflation, and check 
whether this is consistent with the equilibrium solution of the (modified) Boltz- 
mann equation. As was already mentioned in the introduction, the possibility 
of a particle- or string-production-driven de Sitter phase has recently attracted 
some interest ([|l9[, @, [^, 0, |]l2|). We shall show here that this feature is 
not consistent with the kinetic theory of a Maxwell-Boltzmann gas in collisional 
equilibrium. 

It is obvious from ( [T7| ) and (^) that the condition H = Hq = const for 
exponential inflation is 

7r = -{p + p) . (64) 



On the other hand, by ( pig ) an d (p^ ), the effective viscous pressure is given in 
terms of u and r. Using either ([42D or ( |36|) for u and combining (|l^ ) and (p^) 
fixes T to be 



_i / z nm \ ^ 



' n n 



1 dp 



. (65) 
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Solving p5| ) with (pl|) and (22) for r ^ provides us with 



which differs from (47) by the nonvanishing s-term. Ehminating r ^ by (|6^) 
yields 

' = -3d-p^ 

for the change of the entropy per particle. The condition for exponential in- 
flation implies a negative s I As we shall show below, this decrease in the 
entropy per particle is compatible with the second law of thermodynamics only 
up to time scales of the order of the Hubble time. For larger times exponential 
inflation will turn out to be forbidden thermodynamically. 
Instead of (^) we now have 

Sf^ = nsT + ns (68) 

for the entropy production rate. The particle production rate F itself depends 
on s according to 

l--^)e-s, (69) 



3 dp, 

as follows from (|2l] ) with (^6|) and (^). A negative s enlarges the production 
rate F compared with the case s = 0. Inserting (|67|) into (|69|) it follows 



i + ^^^)e. (70) 



3 dp p 

The entropy production density becomes 

nfi p dp 

^;. = -0-^^5^0' (71) 

where the expression (|39| ) has been used, is the chemical potential with 
p/T = a. 

Equations (|l3|) and (|70|) imply 

^ = ^l^^e , (72) 

n 6 dp p 

i.e., the number density increases during the period of exponential inflation. 
This is a natural consequence of the fact that the particle production rate ([70|) 
is larger than the expansion rate. 

We point out that the considerations of this section are not restricted to 
massive particles. Different from the adiabatic case of the previous section 
that did not allow for the production of ultrarelativistic particles with z — > 0, 
there is no corresponding restriction in an exponentially inflating universe. The 



reason for this is a simple one: Since (35) does not depend on r, any choice 



of the latter may be compatible with the equilibrium conditions. It was the 
adiabaticity condition of the previous section that introduced a proportionality 
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between and u. Only by virtue of this additional requirement r was related 
to ly. Since v neccessarily vanishes for z ^ the condition for adiabatic particle 
production implies — > as well. In the present case of a constant Hubble 
rate the corresponding condition on r is independent of v. Consequently, even 
for vanishing a nonzero particle production due to r^^ > 0, that does not 
affect the behaviour ( p5| ) of the temperature, is possible. 

With pr = Pr/3, where the subscript 'r' refers to radiation, equation ([7^) 
yields 



— — -Ho 
rir 



(73) 



equivalent to rir oc R. Since according to ( p6|) we have Tr oc R^^ always, this 
behaviour is consistent with pr = SrirTr = const in the de Sitter phase. 

With i/j. = for radiation, the entropy per particle changes according to 
(H) with (H) and (H) as 

Sr = -- . (74) 



and 



Taking into account (p9|) and Pr+Pr = 'irirTr, the expression (74) is compatible 
with ar = Pr/Tr and the first equation (p4[). Furthermore, in the ultrarelativis- 
tic case the relations 

(75) 
(76) 

and (^) , the condition 
(77) 



5fl 



Tr 



rir , 
\Sr 



1) 



are valid. With Vr = for pr = pr/3 and using (|T 
(|6^) for exponential inflation reduces to 

Tr = -Hn^ 



The characteristic time scale for particle production is a quarter of the Hubble 
time. Combining (|39| ) and (|77|), the entropy production density (^) may be 
written as 

Pr 



S^^ = AurHo ( 3 



The quantity Pr/Tr is given by |2C] 



Pr_ 

Tr 



In 



vr 



(78) 



(79) 



While Pr/Tr is constant for F = 0, this is no longer true in the present case. 
From dl), (|7|) and (|^) we have 



^ =4Fn 



and, after integration. 



Pr_ 
Tr 



Pr_ 

Tr/0 



+ AHo {t - to) 



(80) 



(81) 
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The first term on the r.h.s. corresponds to the initial value at t = to- The 
ratio fir/Tr grows linearly with t. From (|8l| ) and ([78| ) it is obvious that the 
fulfillment of the condition 5*.° > is only guaranteed for 



t-to< -Hq ^ 



o V-tr/ 0. 



(82) 



After a time of the order of the expansion time the second law becomes violated, 
i.e., the corresponding process is forbidden. 

Denoting the opposite limiting case z ^ 1 for nonrelativistic matter by the 
subscript 'm', we find from (72) the expression 



rirr, 2 m 



(83) 



for the time dependence of the particle number density, i.e., in the limit z ^ 1 
the number density is almost constant. The change in the entropy per particle 
is 



e 

2" 



The quantities and Um are given by 



z 1 
3 2 



e 



and 



M 



e 



(84) 



(85) 



(86) 



respectively. Although » and | u \ M/n-m ^ 0, the terms proportional 
to z cancel in the expression (pT|) for T, yielding 



0, 



(87) 



in agreement with the corresponding result from ([70[). Introducing pm = "11^111+ 
SrimTm/'i andpm = nmTm into (^), the expressions (84) and (|85|) are consistent 
with (fim/Tm) = ctm = T""^, as required by (|23). The growth rate of Om is 
much larger than the particle production rate (|87|). 

The entropy production density in the case z ^ 1 becomes 



na 



5 

+ 2 



The chemical potential, the number density and z are related by 

27r2 



exp [fi/T] 



-n . 



m?TK2 {z) 

Using for K2 {z) the asymptotic representation for z S> 1 



K2 {z) ~ exp [-z] , 



(89) 



(90) 
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one finds 



7^ = — in 



mT„ 



TLr, 



(91) 



For r = one lias rim oc R and oc R and the quantity z — /j^m/Tm is a 



constant. In the present case we find from (39) and (84) 



Integration yields 



After a time interval 



2 " 



^ --^o(t 

7m, / 2 



to) 



-1 



1 + 



(92) 



(93) 



(94) 



the condition 5" > becomes violated, i.e., the production of massive particles 
is not consistent either with exponential inflation. 

These results show explicitly, that the production of Maxwell-Boltzmann 
particles in equilibrium is compatible with a de Sitter phase only on time scales 
of the order of the Hubble time. 

In order to avoid misunderstandings we point out that this statement does 
not necessarily imply that all the above mentioned fluid cosmological scenarios 
([H]) 0) Hi! 0' EH) wrong. We have only excluded the case that the 
condition ( |64D for exponential inflation and the equilibrium condition ( p6| ) are 

are dropped, 



fulfllled simultaneously. If the equilibrium conditions (24) - ( psj) 
the conditions for exponential inflation are less restrictive. Retaining the equi- 
librium conditions corresponds to the so-called 'global' equilibrium, while drop- 
ping them and having nevertheless an equilibrium structure of the distribution 
function corresponds to the 'local' equilibrium case (|0], ||l5|). Taking into 
account deviations from either of the equilibrium distribution functions implies 
an imperfect fluid universe which requires separate investigations (|2^, |28|). 

Instabilities of the de Sitter universe due to quantum processes ( p^, ^ ]) or 
due to viscous pressures, modelling particle production phenomena ([^]) are 
well known features. A breakdown of the de Sitter phase due to a violation of 
the condition S'.'J^ > as in the present paper may be regarded as thermody- 
namical instability. If the classical evolution of our model universe is assumed 
to start in a 'global' equilibrium de Sitter phase, our results imply the latter 
to become thermodynamically unstable after a time of the order of the Hubble 
time. On the basis of our setting there exist two qualitative possibilities for 
a subsequent evolution of the universe in accordance with the second law of 
thermodynamics. According to the flrst possibility the universe may leave the 
'global' equilibrium state and stay in a 'local' or nonequilibrium de Sitter phase 
([p7|, |2^). Alternatively, the universe may leave the de Sitter phase altogether 
but remain in 'global' equilibrium, being characterized, e.g., by adiabatic par- 
ticle production, including the possibility of power-law inflation as in section 
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3 above. Whether one of these scenarios will actually be realized, can only be 
decided once reliable knowledge about the functions v and r, that characterize 
the particle production, is available from the underlying quantum physics. 

5 Conclusions 

Within an 'effective rate model' we have studied the equilibrium conditions for 
a Maxwell-Boltzmann gas with variable particle number. Our main findings can 
be summarized as: 1. Radiation and nonrelativistic matter may be in collisional 
equilibrium at the same temperature provided, the number of matter particles 
increases at a specific rate. 2. Collisional equilibrium for massive particles 
with adiabatically increasing particle number is only possible in a power-law 
inflationary universe. 3. Exponential inflation is consistent with the 'global' 
equilibrium conditions for a Maxwell-Boltzmann gas only for a time interval of 
the order of the Huble time. For larger times a 'global' equilibrium de Sitter 
phase violates the second law of thermodynamics. 
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